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A  METHOD  OF  DELETING  THE  OUTLIERS  IN  OPTICAL 
OBSERVATIONS  OF  ARTIFICIAL  SATELLITES 

Wu  Lianda  Jia  Peizhang 


ABSTRACT 

...  ar^-^.c^e  discusses  in  detail  deletion  problems  associated 

with  outliers  in  optical  observations  of  artificial  satellites .  We 
first  of  all  take  complex  nonlinear  problems,  making  use  of  initial 
orbits  and  turn  them  into  linear  models  convenient  for  digital 
processing.  Use  is  made  of  two  quantities— equations  of  time  at  as 
weli  as  orbital  plane  normal  deviation  40  in  order  to  analyze  and 
make  determinations.  In  conjunction  with  this,  a  type  of  L- 
estimation  target  function  opting  for  the  use  of  trimmed  least 
squares  is  put  forward  to  carry  out  methods  associated  with 
deletion  of  outliers.  The  breakdown  point  associated  with  these 
methods  is  50%.  At  the  same  time,  computational  examples  are 
given.  In  conjunction  with  this,  discussions  are  carried  out  with 
regard  to  initial  orbit  accuracies  and  deletion  criteria. 
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In  this,  R  is  survey  station  geocentric  range: 

R  is  survey  station  geocentric  vector; 

Z  =  90*-h  is  satellite  zenith  range; 

1  is  satellite  obervation  direction,  and 

/  -  J?(z,180°  -  S)J?(y, 90°  -  (p)f„ 

/0  —  (  cos  h  cos  A ,  —  cos  h  sin  A ,  sin  A)* , 


Here,  S  is  survey  station  local  sidereal  time; 

<f>  is  survey  station  latitude. 

Use  r  to  carry  out  coordinate  transformation.  Take 
coordinates  and  transform  them  into  a  right  hand  coordinate  system 
which  is  composed  of  xyz  axses  with  an  x  axis  pointing  toward  the 
orbital  ascending  node,  and  a  z  axis  pointing  in  a  direction  normal 
to  the  orbital  plane .  One  then  has : 

rN  “  R(.x,i)R(z,Q)r,  (3) 


R(x,6 ) ,  R(y,0),  R(z,0)  in  equation  (2)  and  equation  (3)  are 
coordinate  rotation  transformation  matrices  (expressions  omited) . 
i  and  Q  are  satellite  angle  of  inclination  and  ascending  node 
longitude . 

Assuming  ,  _  , 

rN  «  U,y,z)r,  (4) 

«o  —  <g_l(y/*), 

one  then  has 


Here,  uo  is  the  satellite  longitude  angle  corresponding  to  r.  This 
is  different  from  u  corresponding  to  time  t  radicals  E.  On  the 
basis  of  binary  problem  formulae,  it  is  not  difficult  for  us  to 
calculate  the  satellite  level  longitude  A  o  .  Then, 


In  this,  .  ... 

A  =  M  +  w  is  satellite  level  longitude  associated  with  time 

n  is  satellite  level  motion  n  =  and  m  =  GM. 

Besides  this, 


Ad  —  sin -1(z/r) 


(6) 


AT  and  A0  defined  in  equations  (5)  and  (6)  are  different  from 
the  observation  quantities  A  and  h.  They  already  correspond  to 
satellite  geocentric  movements.  Speaking  in  terms  of  observations 


associated  with  a  station  ring  (short  arc),  so  long  as  initial 
orbital  radicals  are  not  very  wrong,  they  can  then  be  expressed 
with  the  use  of  linear  models: 

AT,-  —  a#  +  <*i0  i  —  *o)  +  Ci  > 

A 0;  —  bt  +  tiOi  —  O  +  ’ll  >  0  “  1  >2  t . ^0  ’  / 


In  this,  £i  and  r|i  correspond  to  observation  errors.  We  assume 
that  they  are  zero  mean  values — mutually  independent  random  values 

associated  with  normal  distributions. 

If  we  are  able  to  precisely  estimate  ao,  al.  bo,  and  bl,  it  is 
then  possible  to  make  use  of  the  magnitude  of  I AT,  -  ao  -  al(tl  - 
to)|'|A0|  —  bo  -  bl  (tl  -  to)|  as  well  as  the  magnitudes  of  al,  bo, 


and  bl  themselves  in  order  to  carry  out  outlier  deletion.  As  a 
result,  problems  associated  with  outlier  identification  and 
deletion  can  be  summarized  as:  already  knowing  measurement  values 

{ATi}  and  {A0i}  (i-1,2 - N),  as  far  as  certain  outlier  values 

contained  iithin  them  are  concerned,  how  to  solve  for  robust 
estimations  of  ao,  al,  bo,  and  bl. 


II.  ROBUST  ESTIMATION  METHODS 


Assume  that  the  already  known  guantities  {yi}, 
obey  the  linear  model  described  below: 

y.  -  +  btOi  -  h)  +  e‘> 


(i  =  1,2, ...N) 


(8) 


in  the  equation,  {ei>  (i  -  1,2 . N)  are  observation  errors.  With 

reoard  to  good  observation  quantities,  (ei)  can  be  assumed  to  be  a 
re?S  mean  value.  Mutually  independent  random  quantities  associated 
with  normal  distributions  have  variances  of  ou  ,  When  {yi}  has 
certain  outliers  in  it,  the  robust  estimation  methods  associated 
with  ao  and  bo  are  very  numerous  [1-5].  We  opt  for  the  use  °f  L 
estimation  methods  associated  with  trimmed  least  squares.  The 
basic  idea  is  to  make  target  function: 


u 


Nfi+l 

S  (rj 

<»  1 


(9) 


mm, 


in  this  NB  is  a  whole  number  (1/2^B<D-  (rNl)  represents  residual 
-irmences.  It  is  possible  to  demonstrate  that  the  estimates 
in  question  are ’ consistent  estimations.  The  breakdown  point  is  e* 
i  1  -  B [ 3  ]  •  If  we  select  B  =  1/2/  then  the  breakdown  point  is  50%. 

If  we  are  able  to  sample  and  obtain  a  pair  of  good  observation 
quanitites  yi,  yj,  then: 


an  f  (y,  +  y.)/2» 


(10) 


This  is  nothing  else  than  one  type  of  relatively  reasonable 
estimation  of  ao  and  bo.  Because  of  this,  the  problem  is 
summarized  as:  how  to  carry  out  sampling  to  make  it  possible  to  at 
least  be  able  to  obtain  a  set  of  good  observation  quantities.  The 
clumsiest  method  is  to  take  N  pieces  of  data  two  by  two  in  pairs. 

It  is  possible  to  obtain  individual  samples.  It  is  then 

possible  to  obtain  individual  target  functions.  If  there 

are  M  pieces  of  good  data  in  N  pieces  of  data  (M  >  N/2),  then,  in 

C\,  individual  target  functions,  we  can  have  Ci 

individual  relatively  good  estimations.  Making  target  functions  U 
adopt  minima  is  nothing  else  than  robust  estimations  of  ao  and  bo 
associated  with  our  final  solution. 

In  order  to  save  on  calculations,  we  certainly  do  not  need  to 
calculate  r 1  individual  target  functions.  Our  method  is: 

Let 

rl,  10, 

'  m  l/NT(tf/5),  JV>10, 

K-  INT(Nlm), 


Taking  each  individual  datum  and  each  interval  m  individual  datum 
which  follows  it  and  mating  them  up  two  by  two,  the  number  of 
iterations  of  calculations  associated  with  this  type  of  target 
function  is : 

K*[N-m(K+l)/2] 

iterations.  When  N  is  relatively  large,  it  is  much  smaller  than 
£*  .  of  course,  this  type  of  method  is  still  not  the  most 

economical  method.  However,  with  regard  to  when  N^50,  the  method 
is  realistically  feasible.  Moreover,  it  is  possible  to  guarantee 
obtaining  quite  a  large  number  of  sets  of  accurate  data  samples. 
Thus,  it  is  possible  to  select  from  among  them  relatively  good 
estimations  associated  with  ao  and  bo. 

In  the  ordinary  course  of  events,  ao  and  bo  which  are  obtained 
in  this  way  are  also  capable  of  being  further  refined  using  M 
estimations.  However,  due  to  the  fact  that  our  purpose  is  not  to 
find  the  most  precise  ao  and  bo  but  is  to  delete  outliers,  as  a 
result,  we  did  no  further  refining.  Practical  realization 
demonstrated  that  ao  and  bo  which  are  solved  for  in  this  way  are 
already  quite  effective  with  respect  to  deleting  outliers. 


III.  DELETION  OF  OUTLIERS 


{ATi},  {A0i}  opt  for  the  use  of  linear  models  which  are  shown 
in  equation  (7).  Using  the  methods  described  in  sections  above,  we 
are  then  able  to  obtain  robust  estimates  for  ao,  al,  bo,  bl  /116 
Because  of  this,  methods  associated  with  the  deletion  of  outliers 
can  be  summarized  as: 

1.  Random  Outlier  Deletion 
Calculation  residual  errors: 


m  —  AT;  —  [a#  +  «,(*,•  —  f,)] , 

re 1  -  A 9i  -  [*.  +  *,(*,•  -03,  *-1,2,; . N, 


delete 

\rn\  >  0.5s, 

.  I  r0,  |  >  0.03° 

outlier  values. 

2.  Deletion  of  System  Outlier  Values 

System  errors  given  rise  to  by  timing  or  instrument  reading 
drift  or  jumps  as  well  as  errors  given  rise  to  by  observations  of 
erroneous  targets  will  in  all  cases  be  reflected  in  ao,  al,  bo,  bl. 

ao  reflects  time  system  differences  as  well  as  errors  along 
tracks  given  rise  to  because  of  initial  orbital 
inaccuracies ; 

al  reflects  satellite  speed  errors,  including  errors  given 
rise  to  by  inaccuracies  in  satellite  periods, 
eccentricities  e,  and  perigee  angles  w; 

bo  reflects  errors  associated  with  orbital  ascending  node 
longitude  R  as  well  as  errors  given  rise  to  because  of 
surface  height  imprecision  (corresponding  to  satellite 
geocentric  range  r); 

bl  reflects  angle  of  inclination  errors  as  well  as  errors 
given  rise  to  by  r  (corresponding  to  a,  e,  and  w) 
inaccuracies . 

Due  to  initial  orbital  errors,  influences  on  ao  are  very 
large.  We  cannot  use  ao  obtained  by  single  station  ring  data  in 
order  to  precisely  identify  time  system  errors.  As  a  result,  we 
did  not  use  ao  in  order  to  delete  system  outliers. 

Due  to  the  fact  that  our  initial  orbits  had  a  certain 
precision,  as  a  result,  we  are  able  to  use  al,  bo,  and  bl  in  order 
to  delete  system  outliers.  The  method  is  as  follows. 

Delete  I*, I  >0.3° 

|i,j  >  0.1° /minute 
|  ax  |  >  2  s  /minute 

station  rings . 

These  criteria  are  given  on  the  basis  of  initial  orbit 

£ 


precisions  which  we  are  capable  of  obtaining  as  well  as  through 
experimental  test  calculations  of  orbits.  Below,  we  will  discuss 
these  further. 

IV.  CALCULATION  EXAMPLES 

The  radicals  for  a  certain  satellite  on  27  September  1990  0?0 

are: 

a  *  1.09776422  (RE)  e  =  0.004334 

i  =  82 . °  584  l  =  103 . *507 

Q  =  280.722  A  =  68. “672 

P  =  -0.0004  minutes/celestial  ring 

As  far  as  a  certain  survey  station  on  26  and  27  September  is 
concerned,  a  ring  of  various  data  obtained  is  as  follows.  /117 

26  September 


T 

A 

h 

h 

m 

s 

1 

14 

9 

45.68 

277°  0 ' 

17*29 

2 

14 

9 

56.26 

280® 10 ' 

17*46 

3 

14 

10 

25.48 

288  ®  46 ' 

18*09 

4 

14 

10 

30.12 

290°  09 ' 

18*10 

5 

14 

10 

34.94 

291°  36 ' 

18*09 

6 

14 

10 

39.01 

292®46 ' 

18*08 

7 

14 

10 

44.27 

294®22 ' 

18*05 

8 

14 

10 

49.84 

296*00' 

18*01 

9 

14 

11 

40.51 

310*18' 

16*30 

10 

14 

11 

44.93 

311*27' 

16*18 

27  September 

T 

A 

h 

h 

m 

s 

1 

22 

57 

46.12 

35*21* 

17*59 

2 

22 

57 

53.10 

37*07' 

18*36 

3 

22 

58 

28.12 

46*52* 

21*19 

4 

22 

58 

35.85 

49*19' 

21*50 

5 

22 

58 

45.36 

52*24' 

22*26 

6 

22 

59 

03.10 

58*28* 

23*20 

7 

22 

59 

19.21 

64*16' 

23*56 

8 

22 

59 

38.59 

71*24' 

24*14 

9 

22 

59 

53.23 

76*52' 

24*09 

10 

23 

00 

10.95 

83*22' 

23*44 

11 

23 

00 

20.51 

86*47' 

23*23 

12 

23 

00 

37.19 

92*32' 

22*32 

13 

23 

00 

58.96 

99*25' 

21*06 

14 

23 

01 

21.18 

105*46' 

19*24 

15 

23 

01 

46.56 

112*05' 

17*18 

16 

23 

02 

11.23 

117*22' 

15*12 

ic 


During  test  calculations,  AT  and  A0  sequences  obtained  from 
calculations  associated  with  26  September  data  points  1,2,3, 4,  time 
additions  of  0.7s  and  27  September  data  points  1,2,5,7,15  A  plus 
20'  and  points  3,9  A  minus  20 'are  as  follows: 


1 

AT(min) 

-0.0009 

26  September 

* 

A0 ( deg . ) 
0.0122 

2 

-0.0026 

* 

0.0113 

3 

-0.0058 

* 

0.0015 

4 

-0.0061 

* 

-0.0043 

5 

0.0057 

-0.0002 

6 

0.0028 

-0.0032 

7 

0.0045 

-0.0019 

8 

0.0031 

-0.0036 

9 

-0.0014 

-0.0086 

10 

-0.0023 

-0.0067 

1 

AT(min) 

0.0083 

* 

27  September 

A0 ( deg . ) 
0.0500 

2 

0.0114 

* 

0.0430 

3 

-0.0189 

* 

-0.0312 

4 

-0.0028 

0.0026 

5 

0.0127 

* 

0.0160 

6 

-0.0017 

-0.0031 

7 

0.0140 

* 

-0.0036 

8 

-0.0026 

-0.0085 

9 

-0.0169 

* 

-0.0016 

10 

-0.0017 

-0.0011 

11 

-0.0017 

-0.0055 

12 

0.0003 

-0.0041 

13 

-0.0002 

-0.0007 

14 

0.0029 

0.0052 

15 

0.0044 

0.0069 

16 

0.0214 

* 

-0.0437 

Explanation:  AT  and  A0  sequences  which  have  *  after  them  were 
deleted  during  calculations. 


In  the  calculation  examples  above,  error  in  data  observations 

is  o~ 

2’.  Moreover,  artificially  added  errors  are  approximately  20' 
(~10a).  Calculation  results  clearly  show  that  the  methods  in 
question  are  all  capable  of  correctly  deleting  them.  As  a  result, 
methods  are  successful.  Moreover,  we  also  have  no  difficulty  in 
seeing  that  the  breakdown  points  associated  with  the  methods  in 
question  are  already  approaching  50%  (First  ring  error  data 
accounted  for  40%.  Second  ring  error  data  accounted  for  43.75%.). 
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V.  DISCUSSION 


As  far  as  outlier  deletion  in  artificial  satellite  observation 
angles  is  concerned,  the  problems  which  it  involves  are  very 
numerous — for  example,  the  selection  of  robust  estimation  methods 
how  much  initial  information  is  required  by  methods,  the  level  of 
breakdown  points,  and  so  on.  However,  the  most  basic  problems  are: 

1.  requirements  with  regard  to  initial  orbit  precision; 

2.  rational  selection  of  deletion  criteria. 

These  two  problems  are  closely  related.  Large  amounts  of  test 
calculations  clearly  show  that: 

1.  the  smaller  given  criteria  are — that  is,  our  requiring  the 
deletion  of  relatively  small  outlier  values — the  more  accurate  are 
the  initial  orbital  precisions  required  by  methods; 

2.  when  initial  orbital  precisions  are  relatively  bad,  it  is 
easy  to  give  rise  to  phenomena  associated  with  the  irrational 
deletion  of  entire  rings,  that  is,  deletion  of  al,  bo,  and  bl  is 
not 

due  to  observation  system  error,  but  is  given  rise  to  because  of 
intial  orbit  errors. 

The  appearance  of  these  situations  can  be  understood.  We  are 
only  able  to  give  deletion  criteria  on  the  basis  of  the  actual 
status  of  initial  orbital  precision.  In  situations  associated  with 
transit  theodolite  observations,  errors  in  observations  are 
approximately  2 ' -3 ' .  Due  to  the  fact  that  data  used  because  of 
satellite  intermitence  periods  as  well  as  orbital  measurements  is 
also  very  limited,  e  and  A  determination  errors  associated  with 
satellites  are  relatively  large — in  particular,  when  periods  of 
visibility  have  just  begun.  As  a  result,  when  we  select  criteria, 
it  is  believed  that  e  has  Ae~0.002  and  <■>  has  errors  of  30*  or  more. 

Test  calculations  clearly  show  that,  for  the  criteria  given  in  this 
article,  it  is  possible  to  permit  initial  orbits  to  have  errors 
this  large.  Moreover,  when  there  are  only  ascending  phase  (or 
descending  phase)  observations,  initial  orbit  errors  can  also  be 
somewhat  larger — for  example,  w  can  have  errors  of  around  90°  /119 
.  When  initial  orbit  accuracies  are  inadequate — for  instance,  when 
targets  have  just  been  acquired — obviously,  it  is  only  possible  to 
initially  not  carry  out  outlier  deletion. 

For  the  sake  of  appropriately  expanding  the  range  of 
applications  of  methods — during  the  selection  of  criteria — we 
intentionally  make  suitable  relaxations  of  deletion  criteria. 
Besides  this,  when  observations  show  the  appearance  of  plate  jump 
(for  example,  10  individual  pieces  of  data  have  6  pieces  of  data 
that  show  the  appearance  of  system  errors  and  the  number  of  unusual 
values  exceeds  50%  breakdown  points)  the  occurence  of  erroneous 
deletions  will  also  happen.  The  appearance  of  these  situations 
makes  it  so  that  there  are  still  erroneous  data  (errors  relatively 
small)  which  have  not  been  deleted  during  orbital  improvements 
after  preprocessing.  Among  data  which  has  been  deleted  out,  there 
will  also  be  data  which  was  in  fact  correct.  However,  in  the  final 


analysis,  the  proportion  of  unusual  values  is  relatively  small.  As 
far  as  their  processing  methods  are  concerned,  we  will  study  them 
in  detail  in  another  article. 
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